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A method to investigate acoustic Hawking radiation is proposed, where entanglement entropy 
and mutual information are measured from the fluctuations of the number of particles. The rate 
of entropy radiated per one-dimensional (ID) channel is given by 5* = k/12, where k is the sound 
acceleration on the sonic horizon. This entropy production is accompanied by a corresponding 
formation of mutual information to ensure the overall conservation of information. The predictions 
are confirmed using an ab initio analytical approach in transonic flows of ID degenerate ideal Fermi 
fluids. 
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Hawking's prediction of black hole evaporation [l[ is re- 
garded as one of the most remarkable discoveries of mod- 
ern theoretical physics. Black holes radiate like black- 
bodies with a temperature proportional to their surface 
gravity k. As a consequence of Hawking radiation, black 
holes lose mass and eventually disappear. Hawking's pre- 
diction confirms Bekenstein's conjecture that black holes 
possess intrinsic entropy [3], which would be exhausted 
by radiation during the evaporation process. Black hole 
entropy has inspired the development of new concepts, 
such as entanglement entropy |3|, |4|, which is today of 
increasing relevance for quantum information and renor- 
malization group techniques (5|]. 

Hawking radiation itself does not require gravity to ex- 
ist. It occurs in various analogue models of black holes 
0, in particular, in acoustic black holes, the first sug- 
gested model systems consisting ideally of irrotational 
Euler fluids turning supersonic [T'l. These are also called 
dumb holes as sound generated inside the supersonic 
region (the analogue of the black hole interior) does 
not propagate back into the subsonic region. Based on 
the correspondence between the equations of motion for 
sound near his sonic horizon and those of a quantum field 
propagating near the event horizon of a black hole, Un- 
ruh predicted j?} the emission of thermal sound from the 
sonic horizon. The temperature of the radiated sound is 
given by Hawking's celebrated formula [l[ 
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Hereinafter, the Boltzmann constant k is attached to any 
temperature symbol, while entropy is expressed in units 
of it. The black hole's surface gravity k is replaced by 
the gradient of the velocity of the sound propagating up- 
stream [see Fig. l(a)l, evaluated on the event horizon 



where v is the fiow velocity, which is assumed positive, 
and where c is the speed of sound in the comoving frame. 
A way to visualize ^ is to think at the rapidity of which 
sound wave packets diverge exponentially from the sonic 
horizon, i.e., a local Lyapunov exponent in the geometric- 
acoustics limit (see Fig. 1). 

In this Letter, we focus on entanglement entropy 
production of stationary one-dimensional (ID) acoustic 
black holes. By assuming Hawking radiation (see below). 
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FIG. 1. (a) A fluid turning supersonic is a black hole analogue. 
Solid and dashed arrows represent local velocity vectors in a 
space-time diagram of the upstream and downstream modes, 
respectively. The upstream mode has a special role in the the- 
ory as it has a bifurcation on the event horizon characterized 
by the local exponent k. (b) Hawking particles are created 
on the event horizon in correlated pairs, which populate the 
upstream mode. The evolution of such pairs (e.g. the gray 
and black curves) almost simultaneously generated at sym- 
metric displacements from the event horizon in the far past 
time defines a correspondence between pairs of segments, let 
us say A and B, where mutual information is considered. 
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the rate of which entropy flows from the sonic horizon 
into subsonic and supersonic regions is given by 
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Remarkably, this formula has a similar flavor to Pesin's 
formula for Kolmogorov-Sinai entropy rate and various 
generalization of it, where Lyapunov exponents are cru- 
cial ingredients [§-ll|. Here, however, there is no overall 
production of entropy, and mutual information between 
the two parts separated by the horizon forms in order to 
compensate for the entropy produced in each part (see 
below). The predictions for the entropy and mutual in- 
formation are confirmed by an ab initio analytical cal- 
culation in transonic flows of ID ideal degenerate Fermi 
fluids (a microscopic model). Finally, mutual informa- 
tion and entropy could be determined in an experimen- 
tal simulation of a dumb hole with ID quantum fluids by 
measuring fluctuations of the number of particles. 

Entropy and Mutual Information. — Result ([3]) follows 
by assuming that the collective modes propagating away 
from the horizon are thermally populated at Hawking 
temperature. Each mode with wave vector q and disper- 
sion relation uj{q) carries the entropy s = a;/(e^ — 1) — 
ln(l — e~^), where x — Huj{q)/kThh UM- The entropy cur- 
rent @ is obtained by integrating the entropy multiplied 
by the group velocity, i.e., s ^ |f • This simplifies to 

/o°° ^ ^ implying an independence on the actual form of 
the dispersion, provided this is gapless. Instead, the en- 
tropy distributed on a large probe segment of length L is 
calculated by assuming a local linear dispersion for the 
thermally populated upstream modes and is given by 
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In fact, Eq. ([3]) has the meaning of an average entropy 
current; i.e. S is given by the product of ([3]) and the time 
At — L/\v — c\ necessary for an upstream sound signal to 
cross it, provided that nAt ^ 1 in order that the thermal 
wavelength associated to Hawking radiation 



= 27r|c — v\/k 
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could be well resolved on the length L. 

The idea that entropy is produced in a pure quantum 
system may be puzzling at first glance. In order to un- 
derstand how this fits with a null global von Neumann 
entropy, one has to take into account correlations be- 
tween the inside and outside of the dumb hole (see Ref. 
[3| for the astrophysical black hole case). In fact, a mu- 
tual information associated with such correlated creation 
events forms on the event horizon at twice the rate ([3]) in 
order to exactly compensate for the entropy produced. 

This argument can be made quantitative by recogniz- 
ing that Hawking excitations are created in correlated 
pairs, which depart from the event horizon in opposite 
directions with an average relative time uncertainty of 
order of 1 /k and with speed v — c proper of the upstream 



sound mode. The evolution of such pairs indicates a cor- 
respondence between the location of a Hawking particle 
and that of his partner within an uncertainty of order 
of the local Hawking thermal length l^. This defines an 
approximate correspondence between pairs of segments, 
let us say A and B located on opposite sides of the event 
horizon [see Fig. 1(b)]. It is reasonable to require that 
the information on A plus B is not lost; i.e^ for a pure 
quantum state the entanglement entropy [s, |3| of the 
combined system Sa+b — provided At 3> 1/k. More- 
over, since the information encoded in sound fluctuations 
of the upstream mode on A and B is circa the same, 
also the corresponding entropies, taken as measure of the 
maximal possible information content, are so. Introduc- 
ing the mutual information as Iab — Sa + Sb — Sa+b 
helps in summarizing: if information is preserved 
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between corresponding segments A and B. 

Microscopic model. — The model fluid is a degenerate 
ID ideal Fermi fluid ^3j]. The fluid flows in the presence 
of a repulsive potential Vextix), which acts as a ID ver- 
sion of a de Laval nozzle [l?| flxing the sonic event hori- 
zon at the position of his maximum Knax (for convenience 
X = 0). The degenerate flux of incoming particles is orig- 
inated from a reservoir with upper single particle energy 
/i). sufficient higher than T^nax in order to have a finite 
flow over the barrier. The potential T4xt is assumed suf- 
flciently smooth so that the semiclassical description for 
the fluid and the consequent hydrodynamic description 
are appropriate [l3|- The hydrodynamic description of 
the sound propagation and the associated acoustic wave 
equation with his gravitational analogy are given in Ref. 
[15| . The potential maximum acts as a bifurcation for 
particles propagating at energy e — Vniax, which is char- 
acterized locally by the exponent 



K = -d'^{Vc^t/'m)/dx'^ 
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in a classical description of particle propagation. Quan- 
tum mechanics resolves such bifurcation by a finite tun- 
neling, which for asymptotically smooth potential barrier 
is always consistent with a thermal smoothening of the 
Fermi distribution at temperature fcTbh — hn/iTT. In 
fact, here fcTbh is related to top of the hill transmis- 
sion and reflection coeflicients by jrjp = 1 — jifp = 
{1 -I- exp[27r(e — Knax)/^-^]} ^, which is a known result 
[iGj . It is simple to verify that the acceleration k = d\v — 
c\jdx of the upstream sound mode, which propagates at 
the speed v{x) - c{x) = sign(a;)A/2[T4iax - yc^tix)\lm is 
fixed by the curvature of the potential hill by (O. 

For consistency with the hydrodynamic description, 
the energy difference /x^ — Vmax that corresponds to twice 
the hydrodynamic value of mc? at the sonic horizon, 
should be much bigger than fcTbh in order that quantum 
corrections remain negligible in the semiclassical approx- 
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imation. This is reminiscent of Hawking's theory of evap- 
orating black holes where the back-reaction of the radia- 
tion to the gravitational metric is neglected and Hawking 
temperature is consistently assumed small compared to 
the black hole mass M(? . 

The thermodynamic entropy and mutual information 
can be extracted from the microscopic model by using 
the entanglement entropy as introduced by Bombelli et 
al. in Q and Holzhey et al. in j^] in the context of as- 
trophysical black holes. Entanglement entropy is defined 
relative to a certain space region as the von Neumann 
entropy S{p) — — Trplogp, where p is the reduced den- 
sity matrix of a pure quantum state, made "impure" by 
confining it to certain space region 0, 01 • For noninter- 
acting fermions, the entanglement entropy S is directly 
related to the moments of the Full Counting Statistics 
0. In particular, when this last is Gaussian then S is 
just proportional to the number- fluctuation square 
via 



S = (ttVS) 
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Therefore, the mutual information between the segments 
A and B is proportional to the correlation of the respec- 
tive number fluctuations Iab ~ —{tt^ /3)SNaSNb, and 
the problem is reduced to the integration 
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of the two-particle correlation Jp^^-* (x', x) — {p{x')p{x)) — 
p{x')p{x), which is computed below. 

The many-body wave function of the fluid is the Slater 
determinant that corresponds to particles of mass m com- 
ing from the left and occupying scattering wave- functions 
ij^k with wave- vector k between and — ^ImpijY? . 
The scattering wave-functions "(/"fc are defined asymptot- 
ically as ipkix) = exp{ikx) + rfe exp(— ifcx) for x large 
negative and as ipkix) = tkexp{ikx) for x large positive. 
Far from the potential region, the particles move in a 
constant potential, which for the sake of simplicity is as- 
sumed zero. The escape velocity from the potential hill 
?;esc = y/'^Vnia.x/m corresponds to the asymptotic veloc- 
ity |c — w| of the upstream mode. 

For a Slater determinant the two-body correlation 
function is related to the off-diagonal part of the 
one-body density correlation function = 
(*t(a;)^(a;')) via Sp^^Hx',x) = p{x)5{x - x') - 
\p''^\x, x')\^ , where p{x) = p''^\x,x). For the original 
zero temperature noninteracting Fermi fluid p^^\x, x') = 

/o" If 4'k{x)ij}k{x'). We shall start with the computa- 
tion of the fluctuations in the subsonic region. We as- 
sume A located far away from the event horizon in the 
asymptotic flat region. Then p^^\x,x') results in 
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F[Tr{x' -x)/lii], 
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where kesc = mvesc/fi and F{y) = y/ sinh(y) is a function 
that cuts off distances larger than Ih = 27rVesc/K- For 
large negative coordinates x' and x, p^^^ results in 
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where A — x' — x, kp — (k^ + fcosc)/2 is the local Fermi 
wave- vector and pA ~ k-p/n is the density in A. After 
integration of (fTTj) in A, the entanglement entropy results 
in 
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Note that the second term on the right side, which can be 
rewritten as l/31n(A^yi) represents the vacuum contribu- 
tion (zero temperature) to the entanglement entropy and 
corresponds to the renormalized entanglement entropy in 
a conformal field theory 4] . In the large La limit and ne- 
glecting vacuum-type corrections the entropy Sa is given 



Sa = ttLa/GIh = KAi/12, 



(13) 



where At — La/vcsc- Thus, Eq. (jlSp is a microscopic 
verification of Eq. @ and indirectly of the entropy pro- 
duction rate ([3]). The calculation of Sb gives analogue 
results and is not reported here. 

In the following, the mutual information Iab is eval- 
uated. For this purpose, we shall calculate p'''^\x,x') 
between the opposite sides of the barrier. For large neg- 
ative X and large positive x' the term in the integra- 
tion for p^^\x,x') proportional to tk goes to zero like 
l/{x' — x) and is thus negligible. Retaining the remain- 
ing term proportional to r^t^ results in p^-^^x^x') ~ 
eiko..(x'+x)^2l^cosh[Tr{x' + x)/lii]. Thus, the two-body 
correlation is negative and given by 

5p^^Hx,x') = -l/4/|cosh2 [irix' +x)/lii] 



Finally, the mutual information Iab 



(14) 

'{tt^/3)6Na6Nb 
is obtained by the integration ([9]) of (|14l) and results in 

Iab = KAt/6 (15) 

as leading term for L/ln 3> 1. Thus, the above ab initio 
calculation confirms the information conservation con- 
strain for long corresponding segments. 

Result (HH) for the two-body correlation is very similar 
to a corresponding result for a weakly interacting Bose 
gas obtained by Balbinot et al. in Ref. [l3l [see Eq. (8) 
of their paper] . This correspondence further confirms the 
universal character of Hawking radiation. 

To summarize so far, the production of entropy and 
mutual information can be extracted, in virtu of Eq. ([8]) , 
from the measurement of the particle fluctuations 6N^ 
and SNg and of the correlations SNaSNb, respectively. 
This suggests a method to test the entanglement entropy 
associated to Hawking radiation, that would also be feasi- 
ble with the state of the art in measuring fluctuations. In 
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fact, two recent experiments [18| have demonstrated the 
possibility to extract the temperature kT from the fluc- 
tuations 6N^ in ultracold Fermi gases using the classical 
fluctuation result S'^N ~ NkT/mc^. In our nonequilib- 
rium quantum dumb hole, the classical equilibrium result 
is replaced by 

N 2mc\c — v\ 

that is equivalent to Eq. (|13l) by Eq. ([5]) modulo nota- 
tions. Equation (jl6p can be also derived for ID quantum 
fluids (see supplemental material 

A flnite reservoir temperature fcTrcs calculation of , 
not shown here, leads to 6'^N/N — kT\^h /2mc\c — v\ + 
kT-ccs /2toc|c-|-w|. In deriving this relationship, it is as- 
sumed that kT-ccs and fcTbh are much smaller than mc^ at 
the sonic horizon. Moreover, the segment length should 
be much larger than both Hawking thermal length Zh and 
reservoir thermal length It = h{c+v) /kT-^cs so that flnite 
size effects due to zero point fluctuations are negligible. 
However, in an experimental realization with ultracold 
atoms, it may be necessary to include such contributions 
that correspond to (fT2|) . 

Quantum fluids. — Mutual information may be ex- 
tracted from experiments on ID quantum fluids in a sim- 
ilar way as for ID ideal Fermi fluids. At thermal equi- 
librium and for ultralow temperatures that the speed of 
sound is still close to his zero temperature value, the en- 
tropy is related to 5N^ by 



5= (7r/3C«) 5N^ 



(17) 



where ^ ~ h/mc is the vacuum correlation length and 
n the fluid density. Equation ([T7|) is satisfied also when 
only one branch of the sound dispersion is thermal as one 
can see by comparing Eqs. ^ and ^6]) . Equation (fT7|) 
generalizes ^ and may have a broader range of validity. 
Some progress may be gained by assuming relationship 
(jl7p for the nonequilibrium situation of our sonic black 
hole. This may be the natural way to relate the mutual 
information Iab with the correlated fluctuations between 
two segments A and B. I shall tentatively write 
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where is considered a normalization factor for the 
fluctuations 5N on each homogeneous segment. By re- 
quiring the conservation of information Sa+b = and 
using Eqs. (flSl) and (fTS)) it follows that the particle num- 
ber fluctuations SNa and SNb should be anticorrelated 
according to 
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This expression is indeed consistent with the prediction 
of 17| and generalizes the result obtained in the micro- 
scopic model. In fact, using Eq. (8) of Ref. [T^ for 



the two-particle correlation function Sp^'^\x' , x) and per- 
forming the integration ([9]), it results in Eq. ([19]) as the 
leading term in At. Thus, the above conjecture ([TO]) to- 
gether with Eq. (1161) may be useful for the experimental 
investigation of the sonic analogue of Hawking radiation 
in a ID quantum fluid (see also the supplemental mate- 
rial [l9| for a possible finite reservoir temperature). 



By introducing the concept of entanglement entropy 
production in transonic (ID) quantum fluids, acoustic 
Hawking radiation can be interpreted as a dissipation 
mechanism peculiar of a sonic horizon. 
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